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We numerically and analytically calculate the properties of the bipolaron in an extended Hubbard 
Holstein model, which has a longer range electron-phonon coupling like the Frohlich model. In 
the strong coupling regime, the effective mass of the bipolaron in the extended model is much 
smaller than the Holstein bipolaron mass. In contrast to the Holstein bipolaron, the bipolaron in 
the extended model has a lower binding energy and remains bound with substantial binding energy 
even in the large-U limit. In comparison with the Holstein model where only a singlet bipolaron is 
bound, in the extended Holstein model a triplet bipolaron can also form a bound state. We discuss 
the possibility of phase separation in the case of finite electron doping. 



PACS: 74.20. Mn, 71.38.+i, 74.25. Kc 

There is growing evidence that electron-phonon cou- 
pling plays an important role in determining exotic prop- 
erties of novel materials such as colossal magnetoresis- 
tance and high- Tc compounds Since electrons 

in these materials are strongly correlated, the interplay 
between an attractive electron phonon interaction and 
Coulomb repulsion may be important in determining 
physics at finite doping. In particular, when the electron- 
phonon interaction is local, as is the case in the Holstein 
model, finite Coulomb repulsion leads to the formation 
of an intra-site bipolaron (3|-||] , with an effective mass of 
the order of the polaron effective mass |^ . 

It has been recently discovered that a longer-range 
electron-phonon interaction leads to a decrease in the 
effective mass of a polaron in the strong-coupling regime 
1^,0. The lower mass can have important consequences, 
because lighter polarons and bipolarons are more likely 
to remain mobile, and less likely to trap on impurities or 
from mutual repulsion. Motivated by this discovery, we 
investigate a simplified version of the Frohlich model in 
the case of two electrons, 
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where cj j, creates an electron of spin s and a| creates 
a phonon on site j. The second term represents the 
coupling of an electron on site j with an ion on site ^, 
where go is the dimensionless electron-phonon coupling 
constant. While in general long range electron-phonon 
coupling //(j) is considered ^,|^, we further simplify this 
model by placing ions in the interstitial sites located be- 
tween Wannier orbitals, as occurs in certain oxides 
shown in Fig. (^). In this case it is natural to investi- 
gate a simplified model, where an electron located on site 



j couples only to its two neighboring ions, i.e. I — j±l/2. 
We describe such coupHng with fj±i/2{j) — 1 and oth- 
erwise, and refer to this model as the extended Holstein 
Hubbard model (EHHM). We can view the EHHM as 
the simplest model with longer range than a single site, 
and use it to explore the qualitative change in physics 
in the simplest possible setting. While it is clear that in 
comparison to the Frohlich model, our simplified EHHM 
lacks long range tails in the electron phonon interaction, 
the physical properties that depend predominantly on the 
short range interaction should be similar. For example, 
calculating the polaron energy of the original Frohlich 
model as defined in Refs. one finds that 94% of the 
total polaron energy comes from the first two sites. 

In the case when //(j) = Sij, the model in Eq. (|^) 
maps onto a Holstein-Hubbard model (HHM) (see also 
Fig.(|l|b)). The last two terms in Eq. (|l|) represent the 
energy of the Einstein oscillator with frequency uj and 
the on-site Coulomb repulsion between two electrons. We 
consider the case where two electrons with opposite spins 
{Sz = 0) couple to dispersionless optical phonons with 
polarization perpendicular to the chain. 
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FIG. 1. Schematic representation of the simplified a) ex- 
tended Holstein and b) Holstein model on a chain. Filled cir- 
cles represent electron Wannier orbitals, open circles represent 
ions. Solid lines indicate overlap integral t between Wannier 
orbitals, dashed lines represent nonzero electron-phonon cou- 
pling. 

In this Letter we use a recently developed, highly ac- 
curate numerical technique WiM , combined with a strong 
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coupling expansion to study the simplified EHHM. Our 
main goal is to calculate physical properties such as the 
binding energy, effective mass, isotope effect, and the 
phase diagram of the EHHM bipolaron and compare 
them to the Holstein bipolaron that has been thoroughly 
studied recently [|j . Even though the two models appear 
very similar, we find profound differences between the 
physical properties of bipolarons within the EHHM and 
the HHM. 

The numerical method that we use creates a systemat- 
ically expandable variational space of phonon excitations 
in the vicinity of the two electrons [^,^ . The variational 
method is defined on an infinite lattice and is not subject 
to finite-size effects. It allows the calculation of physical 
properties at any wavevector k. In the intermediate cou- 
pling regime where it is most accurate, it provides results 
that are variational in the thermodynamic limit and gives 
energies accurate to 14 digits for the polaron case and up 
to 7 digits for the bipolaron case. 

To investigate the strong coupling regime of the 
EHHM, we use a Lang-Firsov unitary transformation 
H = e^He-^, where S = go Y.ji, fiUhjsiai - a]). This 
incorporates the exact distortion and interaction energies 
for static electrons into Hq, and leads to a transformed 
Hamiltonian 

H = Ho + T, (2) 
Ho^uj^ a^.Gj - ujgl ^ fi {i)fi (j)n^nj + "-iT^ji ' 

j ijl- j 



where nj — nji 



n,x andg2^5§Ea/;(0)'-/;(0)//(i)]; 

g ~ go for the EHHM. The second term in Ho gives the 
polaron energy, which in the EHHM case is ep = 2ujgo, 
while for the HHM, Cp = w^q. This term also includes 
the interaction between electrons located on neighboring 
sites, a consequence of the non-local electron-phonon in- 
teraction. As noted by Alexandrov and Kornilovitch 
in the strong coupling regime a Frohlich polaron has a 
much smaller effective mass than a Holstein polaron with 
the same polaron energy £p. The reason for lower mass 
in the Frohlich case (as well as EHHM) is that the effec- 
tive electron-phonon coupling that renormalizes hopping 
^2 = ^ep/uj is smaller (in EHHM, 7 = 1/2) than in the 
case of the HHM with 7=1. In the strong coupling 
EHHM polaron, the phonon is displaced on two sites. It 
is identical on one of these sites in the initial and the final 
state after the electron hop, resulting in a smaller mass 
enhancement from phonon overlap. 

In the anti-adiabatic limit where go — > and — s- 00 
with ujgQ constant, the phonon interaction is instanta- 
neous and our simplified EHHM model maps onto a gen- 
eralized Hubbard model 
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with an effective Hubbard interaction U = U — 4:UjgQ and 
V — —2ujgQ. In the case of two electrons an analytical 
solution can be found. As many as three bound states 
may exist: two singlets and a triplet. In the case when 
[7 = there is always at least one singlet bound state. A 
triplet bound state with an energy E = —2ujgQ — 2t^ /bjg"^ 
exists only when ug"^ > t. 

In the strong coupling limit, T in Eq. (^ may be con- 
sidered as a perturbation. In the case when U < 2ujgo, 



the single site or SO bipolaron, defined as (f)so — cJ|Co| |0) , 
has the lowest energy to zeroth order. In this regime 
the binding energy is A = E^-^ — 2ep = U — 'iujg^, 
where E^-° denotes the SO bipolaron energy and £p is 
the energy of a polaron in zeroth order. In the oppo- 
site regime, when U > 2ujg^, the inter-site or SI bipo- 
laron, 4>si^'^ — "TI^'^ot'^U '^ox'-it)!^'*' lowest 
energy. Its binding energy A = -~2ujgQ does not de- 
pend on U, which also leads to a degeneracy between the 
spin-singlet (S=0) and the spin-triplet (S=l) state. This 
simple analysis predicts that a EHHM bipolaron (EHB) 
remains bound in the strong coupling regime even in the 
limit when [/ — > 00. 

It is worth stressing that in the limit ?7 — > 00, sin- 
glet and triplet bipolarons become degenerate. We can 
therefore predict the existence of a singlet and a triplet 
bipolaron, where at finite U the singlet bipolaron has 
lower energy. It is also obvious that the energy of the 
triplet bipolaron should not depend on U. In contrast 
to these predictions, a triplet Holstein bipolaron (HB) is 
never stable, and furthermore in the limit U ^ 00 no 
bound HB exists Q. 

Next, we focus on the effective mass of the EHB in the 
strong coupling regime. First order perturbation theory 
does not lead to energy corrections for the SO EHB. Sec- 
ond order perturbation theory gives 
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where m*^^ = <PE{k)/dk'^ . Equation (^ is only valid 
in the limit when 1/A = 2t/ep — s- and U <ti Cp. In 
the limit of large g and U — 0, cx exp(2ep/a;), 

which should be compared to the HB effective mass that 
scales as nigQ cx cxp(4ep/ti;) 11|. In the strong cou- 
pling regime the EHB should be much lighter than the 
Holstein bipolaron. There is a particularly interesting 
EHB regime when U = ep. In this case the zero order 

S—O 1 

energies of the (^so and (pg^ ' bipolarons are degener- 
ate. Degenerate first order perturbation theory can be 
applied to the spin-singlet EHB in this case, which leads 
to a substantial decrease in the effective mass 
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m*EHB{U^e,)^^e^'^'-. (5) 

The EHB in this regime consists of a superposition 
of E^nd (pg'^^j and moves through the lattice in a 
crab-hke motion. Its binding energy is A = — ep — 
2V2iexp(-ep/2tj). 

In the U oo hmit we apply the second-order pertur- 
bation theory to the SI bipolaron. We take into account 
processes where one of the electrons within the SI bipo- 
laron jumps to the left (right) and then the other follows. 
This leads to 

m^^5(C/ = (^) = ^e^-/-. (6) 

Strong coupling approach thus predicts nonmonotonous 
dependence of the effective EHB mass as a function of U 
as can be seen from different exponents in Eqs. (§,||j6|). 
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FIG. 2. a) The bipolaron inverse effective mass vs. A at 
Lj = 1. The thin full line and thin dot-dashed line represent 
strong coupling results obtained using Eqs. (|) and (|) respec- 
tively, b) The effective mass vs. U sit u — 1 and A = 1.45. 
The thin full and dot-dashed lines represent strong coupling 
results obtained using Eqs. and respectively. 

We next present numerical results. To achieve suf- 
ficient accuracy, we have used up to 3 x 10^ varia- 
tional states. We use units where the bare hopping con- 
stant is t ^ 1 . The ground state energy of the EHB at 
A = 0.5, Lu = U ^ 1, is E = -5.822621, which is ac- 
curate to the number of digits shown. (For the same 
parameters, U — 0, the Holstein bipolaron energy is E 
= -5.4246528.) The accuracy of our plotted results in 
the thermodynamic limit is well within the line-thickness. 
In Fig. (|^a) we present the inverse effective masses of 
the EHB and the HB at [/ = and of the EHB at 
U ^ €p. Our results for the bipolaron mass are in qual- 
itative agreement with results for the polaron effective 
mass by Alexandrov and Kornilovitch In the weak 
coupling regime we find the EHB slightly heavier than 



the HB, while in the strong coupling regime the oppo- 
site is true. Setting the Coulomb interaction to U = Cp, 
the effective mass becomes even lighter, which is a conse- 
quence of the smaller exponent in Eq. (|^). In the strong 
coupling regime (A > 1), we find good agreement with 
our strong coupling predictions in Eqs. (0,^), depicted 
by thin lines. While the absolute values may differ by up 
to a factor of 4 (in the case of U = Ep), the strong cou- 
pling approach almost perfectly predicts the exponential 
dependence (seen as parallel straight lines in Fig. (|^a)) 
of the effective masses on ep = 2tX. 

To obtain better understanding of the effect of on-site 
Coulomb repulsion on the bipolaron effective mass in the 
strong coupling regime, we present in Fig. (^) effective 
masses of the EHB and HB at fixed coupling strength 
A = 1.45 as a function of U. The most prominent find- 
ing is that the EHB is two orders of magnitude lighter 
than the HB when U — 0. While the effective mass of 
the HB decreases monotonously with U, the EHB effec- 
tive mass reaches a shallow minimum near U — ep as 
predicted by the strong coupling approach. At larger 
U > ep we observe a slight increase in the effective mass. 
In the same regime HB effective mass drops below EHB 
effective mass. This crossing coincides with a substantial 
decrease of the HB binding energy and consequently with 
separating of HB into two separate polarons. Numerical 
results for the EHB agree reasonably well with analytical 
predictions for small U, Eq. (^), and also in the limit of 
large C/, Eq. (|). 
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FIG. 3. a) Binding energies A^^'^) vs. U of the EHHM 
(full lines) and the HHM (dashed lines). Corresponding bind- 
ing energies of the first excited states are indicated with thin 
lines, b) Phase diagram of the EHHM (filled circles) and the 
HHM (open circles) calculated at lj = 1. The vertical line at 
A = Ac represents the stability line of the S = 1 EHB. Text 
in the figure applies only to EHHM phase diagram, c) The 
inverse effective mass and the isotope effect of the EHHM vs. 
X at U — 5. Vertical lines represent stability limits of the 
S = and S = 1 EHB (from left to right). 

To gain an insight into the symmetry of the bound 
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EHB state, we have calculated the binding energy 
A'"^^) = -S'bi''^^ — 2£'po where are the ground state 

and the first excited energy of the EHHM or HHM for 
two electrons with opposite spins, Sz = 0, and Epo is the 
ground state energy of the corresponding model with one 
electron. In Fig. (^) we present binding energies of the 
bipolaron ground and first excited states as a function 
of U. An important difference between the HHM and 
the EHHM is that in the former case a critical Uc exists 
for any coupling strength A when the HB unbinds, while 
the EHB remains bound even in the limit U oo when 
X > Xc — 0.76. At small U excited states of both models 
correspond to bipolaronic singlets, spaced approximately 
uj above the ground state. Singlets can be recognized by 
the fact that their binding energies depend on U. As U 
increases, the excited state of the HB unbinds while the 
excited state of the EHB undergoes a transition from a 
singlet to a triplet state which is also bound. 

By solving A^^'^^ (A, Uc) = we arrive at the phase 
diagram {Uc,X) of the EHHM calculated at fixed lu = 
1, presented in Fig. (^). We indicate three different 
regimes. For small A and large U no bound bipolarons 
exist. With increasing A there is a phase transition into a 
bound singlet bipolaron state. Increasing A even further, 
a triplet bipolaron becomes bound as well at A = Ac. For 
comparison we also include the phase boundary of the 
HHM (open circles). Note that only a singlet bipolaron 
exists in the HHM. 

In Fig. (^) we present a cross section through the 
phase diagram in Fig. (||b) at fixed U ~ 5, and plot m*^^ 
and the isotope effect a = dlnmbi/dlnM vs. A (see also 
discussion of the isotope effect in Ref. [||). The effective 
mass increases by approximately a factor of 2.5 from its 
noninteracting value in the regime where only a spin- 
singlet bipolaron exists (between the two vertical dashed 
lines). The increase of the effective mass is followed by 
an increase in the isotope effect. The binding energy (not 
plotted) reaches a value A ~ —OM at A = Ac = 0.76. 

To conclude, we have shown that a light EHB exists 
even in the strong coupling regime with an effective mass 
that can be a few orders of magnitude smaller than the 
HB effective mass at small U. At finite U = ep a regime 
of extremely light EHB is found where bipolaron effective 
mass scales with the same exponent as the polaron effec- 
tive mass. This mobile bipolaron arises as a superposition 
of a (j^so s-iid 8- 0S1 state and it moves through a lattice 
in a crab-like motion. As found in ref. HB becomes 
very light with increasing U close to the transition into 
two unbound polarons dA, U — Uc- Near this transition, 
its binding energy diminishes substantially and reaches 
A = at the transition point Uc- In contrast, EHB can 
have a small effective mass even in the regime where its 
binding energy is large ( in the strong coupling regime 
A approaches A = — e^). Furthermore, EHB remains 
bound in the limit when ?7 — > oo. As a consequence of a 



longer-range electron-phonon interaction, a bound spin- 
triplet bipolaron exists in the EHHM for A > Ac The dif- 
ference between the binding energies of the spin-singlet 
and the spin-triplet bipolaron is proportional to 1/U. In 
the weak to intermediate coupling regime of the EHHM 
(A < Ac and finite U) S = bipolarons exist with sub- 
stantial binding energy close to A ~ Ac, and an effective 
mass of the order of noninteracting electron mass. 

The existence of a singlet and a triplet EHB state has 
important implications in the case of finite doping. As 
was established previously, there is no phase separation 
in the low-density limit of the HHM despite a substan- 
tially renormalized bandwidth Q . The reason is in part 
that a triplet bipolaron is always unstable. The lack 
of phase separation in the low-density limit and in the 
strong coupling regime has a simple intuitive explana- 
tion: a third particle, added to a bound singlet bipolaron, 
introduces a triplet component to the wavefunction. The 
opposite is true in the strong coupling limit of EHHM 
where singlet and triplet bipolarons coexist. In this case, 
the third added particle simply attaches to the existing 
singlet bipolaron and thus gains in the potential energy. 
We therefore expect that the EHHM phase separates in 
the case of finite doping for A sufficiently large. To sta- 
bilize a system of EHHM bipolarons against phase sepa- 
ration, a long-range Coulomb repulsion should be taken 
into account. This prediction is in agreement with recent 
findings by Alexandrov and Kabanov jl^] that state, that 
there is no phase segregation in the Frohlich model in the 
presence of long-range Coulomb interactions. 
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